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ABSTRACT: We study the XY model on a lattice with fluctuating connectivity. 
The expectation is that at an appropriate critical point such a system corresponds 
to a compactifled boson coupled to 2d quantum gravity. Our simulations focus, 
in particular, on the important topological features of the system. The results 
lend strong support to the two phase structure predicted on the basis of analytical 
calculations. A careful finite size scaling analysis yields estimates for the scaling 
exponents in the low temperature phase. 



Introduction 



In recent years considerable progress has been made in understanding the na- 
ture of some simple systems encorporating two dimensional quantum gravity. These 
systems are thought to correspond to string theories out of the critical dimension. 
Two approaches have been followed which are at first sight quite different. The first 
of these employs techniques borrowed from conformal field theory to calculate the 
spectrum of anomalous dimensions in several simple models [1-2]. In the second 
method a regularisation of the continuum functional integrals is made in terms of 
random triangulations [3-5]. Where it has been possible to solve the models ana- 
lytically the two approaches have been in complete agreement. Furthermore, the 
triangulated models may be studied nonperturbatively by numerical simulation. 
Whilst the analytical work gives only information on the scaling behaviour of mat- 
ter field expectation values, the numerical studies allow us to probe the geometrical 
and fractal structure of the worldsheet. The lattice models may also be used for 
matter field central charges c greater than unity where the continuum methods 
break down. 

The marginal case when c = 1 can be realised by employing a single scalar field. 
If the latter model is compactified by imposing that the field be periodic in some 
interval, no exact solution has been found. However there are strong theoretical 
arguments to suggest that the model may exist in two phases; one in which the 
scalar field is non-propagating and is described by pure (c = 0) gravity and another 
corresponding to a massless, propagating scalar field with unit central charge [6] . 

The lattice model is also interesting from a statistical mechanics point of view. 
The model in fiat space (regular lattice) is known to undergo a Kosterlitz-Thouless 
KT phase transition at some critical coupling to a massless low temperature phase 
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with a set of continuously varying critical exponents. This transition is driven by 
the condensation of topological defects (vortices) corresponding to field configu- 
rations with non-zero winding number. It is interesting to ask how this picture 
is modified by the encorporation of lattice (or metric) fiuctuations. Naively, one 
might imagine that performing an annealed sum over lattices might strongly influ- 
ence the effective forces between vortices and hence alter the character of any phase 
transition. Much of the results we present later will be concerned with elucidating 
the role such topological configurations play on a dynamical lattice. 



Discussion of the regular lattice model 



Consider initially the model on a fixed, regular lattice. The latter is taken to be 
the dual to a regular triangulation - a (p^-graph of toroidal topology and possessing 
N nodes. The partition function defining the model is 

Z = $:e^^(/^) (1) 
with the lattice action S {(3) given by 

S{(3) ^ -(3j2^o^i^i-^j) 0ie0---2TT (2) 

{ij) 

Using techniques drawn from conformal field theory it is possible to show that the 
critical region of this theory can be described in terms of two sets of operators with 
anomalous dimensions {An, Am} (see for example [7]). 



Am = l (27r/3eff ) (4) 



The former operators (An) correspond to spin waves, whilst the latter (A^) are 
to be identified with vortices. The effective coupling constant PqQ is determined 
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by details on the scale of the lattice cutoff, for example the lattice type and bare 
lattice coupling {3. However, for large coupling (3 we expect 

/?eff ~ ^ oo (5) 

The vortex operators are seen to be irrelevant operators at large (3qQ ~ /3 and the 
critical theory can be understood in terms of its spin wave content. In this region 
we expect that the theory is massless and the two point function corresponding to 
spin wave excitations behaves as 



However for small j3 this situation is reversed and the vortex operators are relevant 
and dominate the long distance physics. A phase transition is thus expected be- 
tween a region where spin waves essentially exhaust the physically relevant degrees 
of freedom and one where the ground state is populated by vortices. This transi- 
tion will occur when the the anomalous dimension of the lowest vortex operator 
Ai is just marginal (Ai — 2). This yields an estimate for the critical coupling 
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A careful renormalisation group analysis [8] shows that the phase transition at (5c 
is of an unusual type - the free energy and all its derivatives are continuous at f3c 
whilst the correlation length and susceptibility possess essential singularities as the 
transition is approached from the disordered phase. 

^ - aeT^^ (8) 



^^cei0c-pr (9) 

The exponent v takes the value v = 1/2. Furthermore, the dimension of the lowest 
spinwave operator is then 

Ai = ^, (3 = (3c (10) 
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The presence of such a phase transition can be inferred in another way by 
considering the free energy of such a vortex configuration. At large the cosine 
action may be replaced by an equivalent gaussian action, provided we retain the 
periodicity of the 9 coordinates. A simple vortex configuration would be ^ = mcj) 
in an (r, (p) polar coordinate system. This yields an action for the vortex of 

2- fL 



'S'vortex = -/^effTTn In y-j (11) 

where L is the infrared and a the ultraviolet cutoff. Since there is a choice of 
{L/af positions on a cubic lattice at which to place the centre of such a vortex 
(with n = 1) the contribution to the free energy of such a configuration is 

= (7r/3eff - 2) In (^^^ (12) 

Clearly for large /3gff vortices will be irrelevant, whilst for small coupling they will 
dominate, disordering the ground state and forcing a finite correlation length. 

The presence of such non-trivial field configurations is signaled by a non-zero 
vorticity V defined in the continuum by 

V = — I Vd.ds (13) 
271 Jc 

We have employed the naive transcription of this object onto the lattice 

Nl 

y=Y.Y. (^i+l - ^i) 27r (14) 
loops j 

The outer sum is carried out over all the Ni^ loops of the 0^-graph, whilst the 
inner sum corresponds to computing the net change in angle on traversing a given 
loop. In practice, the total vorticity vanishes as a consequence of the boundary 
conditions, so we measure the modulus of the local vorticity in our simulations. 
Clearly, at large where all the spins are locally correlated, this is likely to be 
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a good approximation to the continuum expression, but one might worry about 
the situation at small coupling. This is, of course, the usual problem of defining 
topology on a lattice. The problem is to construct lattice observables which go 
smoothly over into continuum topological quantities. Typically, the lattice oper- 
ators receive large contributions from fluctuations at order the cutoff which have 
no place in the continuum expectation values. One of the standard methods to 
measure topological charge for lattice QCD is by a so-called cooling method [9], 
whereby the action is locally minimised to remove ultraviolet fluctuations. We have 
tested such a method in this situation and found that for a wide range of couplings 
close to criticality the winding number remains constant under cooling giving us 
confidence that we are not seeing lattice artefacts. 



Fluctuating lattice model 



The simulations described later will concern themselves with a nonperturbative 
study of this model when formulated on a dynamical random lattice. As we have 
described this prescription is equivalent (at least when c < 1) to coupling the spin 
model to gravity. Since the XY model has c = 1 at large the KPZ framework [1] 
marginally applies and we can attempt to derive some predictions for this model 
in the presence of such metric fluctuations. Speciflcally, if we take some operator 
O (x) and integrate it over the surface (to get a reparametrisation invariant object), 
the expectation value will scale with surface area as 

2^ n:r^r^\\ . a1-D 



d^x^O{x)) A^oo (15) 

/ gravity 

The KPZ formula then relates this gravitational scaling dimension D to the anoma- 
lous dimension A of O (x) . 

D=^fKl2, c=l (16) 
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Similarly, an integrated correlation function or susceptibility conjugate to O (x) 
would scale like 

/ [ (fx^O{x) j (Px'y^O{x)\ ~^(2-2^) (17) 

J I gravity 

To adopt a conventional normalisation of the susceptibility we would divide by a 

power of the area to get 



X 



(0)^a(1-2^)=a(^-^) (18) 



which is to be compared with the usual finite size scaling in the absence of gravi- 
tational fluctuations 

X^^* (O) ~ (19) 

For the gravity coupled system the condition for relevance of an operator O 
is now that its gravitational scaling dimension D (O) = 1, which in the case of 
the vortex operators O = Ai yields the condition that Ai = 2. This is precisely 
the same condition that determined the critical coupling in flat space. Thus this 
continuum analysis would predict that the gravitational dressing of the vortex 
operators does not change the position of the phase transition. However it is not 
clear that the transition is still of KT type. We will address this question again 
when interpreting our results. 

The simulation of this dynamical model is effected by means of two forms of 
update. To simulate the sum over spin conflgurations we have employed the cluster 
update pioneered by Wolff [10]. This proves to be a very efficient method for quasi- 
massless theories like the XY model. The sum over graphs is effected by a local 
procedure detailed in a previous paper [11]. We have simulated lattice volumes 
from iV = 100 to iV = 5000, typically using O (lO^) sweeps of the lattice per 
coupling Errors were assessed by the usual rebinning procedure. 
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Results and analysis 



Since the continuum arguments suggest that the transition on the dynamical 
lattice may be rather similar to the regular (flat space) theory, we chose to do some 
moderate simulations of the latter to compare with the dynamical results. Fig.l 
shows a plot of the specific heat C 



as a function of lattice coupling (3. The plot also shows the topological susceptibility 
Xv determined by 



Clearly, both show peaks in the neighbourhood of /3c ~ 1-6 which may be taken 
as a signal for some form of critical behaviour. However, the specific heat rapidly 
saturates with increasing system volume and hence is not a good observable for 
determining the nature of any phase transition. Furthermore, in the case of the 
dynamical system the peak in C occurs at a substantially different coupling than 
our best estimate of (3c. Similarly, the spin susceptibility increases monotonically 
with (5 and yields no information on the position of any phase transition. 

In contrast, the topological susceptibility being determined by the degrees of 
freedom which are believed to drive the phase transition, is a much more reliable 
indicator of criticality in the system. Indeed, we have taken the peak in Xv ^is 
our best estimate of the critical point. From the largest regular lattices we have 
considered we estimate the critical lattice coupling to be (3c = 1-7(1). As we detail 
below this produces estimates for the critical exponents on a regular lattice which 
are in agreement with analytical calculations. 

Fig 2. shows the vorticity V for a 3000 node dynamical lattice over a range of 




(20) 




(21) 
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couplings spanning the critical region. The absence of any discontinuities makes 
it unlikely that there are any first order transitions in the system. However the 
topological susceptibility is much more interesting and is shown in fig. 3 for the 
full range of system sizes we studied. The first point to note is that the critical 
coupling, as dictated by the peak in has substantially shifted from its value 
on the regular (flat space) lattice /3c = 1-7 to a value which we estimate from the 
peak on the largest lattice to be (3c = 2.3(1). Note that this is in agreement with 
an earlier study which utilised an entirely different method to derive a value for f3c 
[12]. 

However, at first sight, this renormalisation of the coupling would appear to be 
at variance with the continuum calculations which indicate that transition coupling 
is not changed by the metric fluctuations. In an effort to cast some light on this we 
simulated an XY model on a single quenched lattice deriving from a simulation of 
pure gravity. In this case the peak in Xv again occurred in the vicinity of (5c = 2.3. 
We interpret this situation in the following way. The transition occurs at some 
value of /3gff which is identical on both flat and fluctuating geometries. However the 
corresponding bare lattice coupling may be different in the regular and dynamical 
lattice cases (since their local lattice structure is radically different). Furthermore, 
it is then conceivable that the KPZ predictions (eqn. 18) for the fluctuating lattice 
will still hold true at the appropriate critical point (and indeed our later results 
will lend support to this conclusion). 

The second observation to be made from fig. 3 is that Xv shows strong finite 
volume behaviour for large (3. This is also true of the model restricted to a regular 
lattice and comes about for the following reason. Consider the regular lattice 
topological susceptibility Xv^^ as a (euclidean) time t integral over a vortex-vortex 
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correlator. 

Xr = EEK-0l«)|V^^* (22) 
t s 

The states s are eigenstates of the lattice Hamiltonian with energies Eg. Clearly 
the integral will be dominated by the lightest state sq which couples to the local 
vortex operator vq. It will be this state which will govern the finite size scaling 
of the susceptibility. If we regard the matrix elements {vq\s) now as expectation 
values in a path integral, we can use general symmetry principles to determine 
which of them, in principle, are nonvanishing. Consider the overlap of a spinwave 
state s = e^^ on the local vorticity vq. If we perform a parity transformation of the 
form 

xi — > —xi, X2 — > X2, 0^9 (23) 

the vorticity changes sign but the spin wave state is invariant. Thus, provided 
parity is not spontaneously broken, the spinwave and vortex sectors are completely 
orthogonal. However if we take as operator the modulus of the local vorticity (which 
is the case in our numerical simulations) then both operators are invariant under 
this parity transformation and hence we would expect on these general grounds 
that the lightest state coupling to such an operator would again be the simple 
spin wave state. Thus the finite size scaling of the topological susceptibility will be 
governed by the same critical exponent Ai as for the spin susceptibility. Similarly in 
the gravitationally dressed case of the dynamical ensemble the two susceptibilities 
will scale with the same exponent. A measurement of both quantities gives us 
independent estimates of this scaling exponent. 

This picture is borne out by our numerical data. Fig. 4 reveals an analysis 
of the finite volume behaviour of the spin susceptibility Xs for both regular and 
fiuctuating lattices. As we emphasised, we expect that its finite size scaling will be 
dominated by the lightest spin wave state (n = 1). 
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The regular lattice data shown derive from runs at the transition point /?c = 1-7 
We fit our data from lattice sizes N = 968 to N = 5000 to the form 

X ~ a + bN"" (24) 

using a simple nonlinear least squares procedure. The fit yields a;^^* = 0.855(3) 
with x^/dof = 0.4 with 55% confidence. This number is to be contrasted with 
the prediction a;^^* = 0.875 dictated by the usual KT and CFT arguments. How- 
ever the fiuctuating lattice gives the fit a; = 0.53(13) with x^/dof = 0.04 at 85% 
confidence. This, in turn, lies remarkably close to the KPZ prediction of a; = 0.5 
derived from eqn. 18 with Ai = 1/8. 

Similar conclusions can be drawn in the case of the topological susceptibility 
Xv (fig- 5). Here, the regular lattice fit gives a;^^* = 0.843(31), x^/dof = 0.02 with 
a confidence level of 89%, which is consistent with the spin susceptibility measure- 
ment and compatible with the theoretical prediction. The dynamical lattice yields 
a; = 0.64(15), x^/dof = 0.1 at 79% confidence. This again is in reasonable agree- 
ment with the KPZ number of 0.5 and is inconsistent (at two standard deviations) 
with the fiat space scaling exponent uj = 0.875. 

Clearly then, we have seen very good agreement with the usual KT assignments 
in the case of fiat space. This gives us confidence that our method for locating 
the transition coupling is a good one. Furthermore, the dynamical exponents are 
dramatically different and lie intriguingly close to the naive KPZ prediction of the 
continuum theory, even the critical lattice coupling for the dynamical ensemble is 
significantly different from its fiat space value. 

As an independent check we have examined the spin correlator g (r) for the 
regular lattice at our best estimate for Pc = 1-7. Fig. 6 shows an 'effective' Ty-plot 
determined from the spin correlator on a lattice with N = 2888 by assuming power 
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law behaviour 



rj (r) 




(25) 



We expect finite size effects to be important wlien r > 15 and lattice spacing errors 
to be present for small r, but the figure clearly shows a rather stable plateau for r] 
within this region. Indeed the plot of fig. 7 allows a fit for rj in this range yielding 
rj = 0.247(1) with x^/dof = 1.4. This is to be compared with the predition of 
r] = 2Ai = 0.25 from KT theory. The 'effective' mass computed in the analogous 
way on the assumption of an exponential behaviour for g (r) is shown also in fig. 
6. Clearly there is no stable region for fitting here and we interpret this as strong 
indication that at the critical point of the regular lattice XY model there are power 
law correlations. 

Contrast these conclusions with fig. 8 which shows the same quantities com- 
puted with the ensemble of random graphs. In this case we have to be more careful 
in defining the correlation function g (r). On the lattice we define g as follows 



The function dij is just the geodesic distance (minimal length walk) on a given 
graph between sites i and j, whilst 



We average over spin configurations — e * and graphs G]\[. For a fixed lat- 
tice this just reduces to the usual definition of the propagator. Notice that the 
continuum version of this quantity is just 




(26) 




(27) 




(fix / d^'du ^\fg' (\) [x] (\) (x') b [d (^g, x, x'^ — r) 



(28) 
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This correlation function is explicitly reparametrisation invariant (if the fields are 
scalars) since the identity and length of the geodesies d (x,x', g) are reparametri- 
sation invariant. Therefore this correlation function is a bona fide observable for 
gravity coupled systems. However, in this dynamical case, the effective power fit 
plot at (3 — 2.3 (fig. 8), is extremely poor, and the exponential fit proves consider- 
ably better. Indeed, the best fit we could produce (fig. 7) corresponded to a two 
exponential fit from distances 3 to 15 of the form 

g (r) = a exp —br + cexp —dr (29) 

with the results a = 0.23(1), b = 0.32(4), c = 0.69(2) and d = 0.048(2) with a 
X^/dof = 0.28 at 97% confidence. The second exponential presumably corresponds 
to an excited state, but at present it is not at all clear why there is this qualitative 
difference in the structure of this correlation function g when the system is coupled 
to gravity. 

We have also attempted to fit the spin susceptibility for the dynamical lattice 
simulations to the KT form of essential singularity. Using the data at N = 3000 
and with Pc = 2.3 as before, we have assumed u = 0.5 and attempted to fit the 
data in the disordered phase. Fig. 9 illustrates this by plotting In Xs as a function 
of {Pc — /3)~^'^- The rounding visible in the plot corresponds to the onset of finite 
size effects, but clearly the data set is not inconsistent with a KT type singularity. 
However, it is perfectly possible to fit the data reasonably well over the same range 
by a simple power fit and objectively it is very difficult to assess which is the 
better. This situation is well known in the case of the ordinary XY model and it is 
only relatively recently that numerical simulations have proven capable of resolving 
this problem [13]. With the lattice sizes employed in this study, such a definitive 
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conclusion is impossible. However we can say that a fit of the form 

X~a(/?c-/?)-^ (30) 

requires, typically, a very large coefficient a ~ 40.0 and power 7 = 3.5 — 4.0. 

To examine the effective gravitational action induced by the coupling of spins 
to the dynamical lattices, we also measured the ffuctuation in the local scalar 
curvature. This is just defined by 

'•'(/') = (E^^) (31) 

where li is the length of loop i and the sum runs over all loops in the (/> -graph. 
Fig 10. shows this as a function of the bare lattice coupling for a variety of lattice 
volumes. Clearly finite size effects are small, and the coupling to matter enhances 
the production of singular geometries. Notice that the the effects of this are almost 
constant throughout the low temperature phase, perhaps suggesting that, as has 
been observed before [14,15], the effective action for gravity may be governed only 
by the central charge; the field content and anomalous dimensions seem to have 
little infiuence on the gravitational sector (at least in those models examined so 
far). The figure also shows the same quantity for two Ising models confirming this 
universality. 

On the same figure we show a quantity W corresponding to the cross correlator 
of the local vorticity (again its modulus strictly) with the ring length. 



W{(3)^^ '--1 (32) 

T.i'{k){\vi\) 

The signal is normalised by the disconnected part to factor off the trivial (3 
dependence of the vorticity. Again a peak is seen close to the initial peak in 
the scalar curvature ffuctuation and significantly to the left of the estimated 
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transition point. This would indicate a significant tendency of vortices to bind to 
curvature defects on the lattice. The latter can be understood qualitatively as a 
consequence of the logarithmic divergence of the free energy of a single vortex. The 
free energy is lowered if the vortex forms on a ring of large length (the ultraviolet 
cutoff a in eqn. 11). This conclusion appears to be essentially unaltered when 
fluctuations are taken into account. 

Finally we have studied the clustering properties of vortices and antivortices. 
To do this we introduce the joint probability distributions P^y (r) and Pvav ir)- If 
we have a vortex at the origin r = 0, these measure the probabilities of finding 
a vortex (in the case of Pyy) or antivortex (Pyav) at some geodesic distance r 
(measured with the respect to the loops). Fig. 11 illustrates this with histograms 
of these distributions in the disordered phase and ordered phase. Clearly the 
vortex-vortex distribution is essentially the same in both cases and corresponds 
to an approximately uniform distribution of vortices over the plane. However 
a dramatic difference is seen in the case of Pvav- In the disordered phase the 
there is a no evidence of any real clustering, whilst in the low temperature phase 
almost all antivortices are bound one lattice spacing away. This is completely 
analogous to the situation in flat space where the disordered phase corresponds to 
a uniform plasma of vortices and antivortices, whilst in the low temperature phase 
the vortices and antivortices exist only in bound pairs. The transition is then 
pictured in this language as the dissociation of these dipole-like pairs, whilst the 
peak in Xv corresponds to a maximal fluctuation between pairs as their effective 
binding decreases to zero. Thus, the qualitative picture of a transition driven by 
condensation of defects persists even when the defects are dressed by gravitational 
fluctuations. 
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Conclusions 



We have presented results for the XY model coupled to a fluctuating lattice 
of fixed topology. These simulations are motivated by the need to understand the 
coupling of 2d gravity to simple matter systems. We have collected data both for 
the standard spin wave observables and also the topological sector. The estimates 
of the critical exponents we derive from finite size scaling analyses are in agreement 
with both flat space KT predictions and the KPZ framework. 

We observe that the fluctuations in the geometry plateau in the low temper- 
ature phase at a value commensurate with a two Ising system which is consistent 
with the assignment of central charge one for this phase. Furthermore we observe 
a strong binding of the vortices to regions of high curvature close to the critical 
point. The distribution functions for vortex-vortex and vortex-antivortex pairs 
conflrm the qualitative picture of the two phases as given by KT theory in flat 
space. The presence of metric fluctuations, whilst dressing operators, seems to 
preserve the nature of the transition. 

These results then conflrm the phase structure argued for on the basis of contin- 
uum calculations, and furnish yet further evidence of the validity of the dynamical 
triangulation prescription (and its computer implementation) in the marginal case 
of c = 1 theories. 

This work was supported, in part, by NSF grant PHY 92-00148 and some of the 
numerical calculations were performed using the Florida State University CRAY 
YMP. We acknowledge fruitful discussions with Ian Drummond and Sumit Das. 
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FIGURE CAPTIONS 



[I] Specific Heat C and Topological Susceptibility Xv for a regular lattice {N = 968) 
as a function of bare coupling j3. 

[2] Vorticity on a 3000 node dynamical lattice vs (3. 

[3] Topological Susceptibility Xv on a dynamical lattice. Lattice volumes are = 100 
(x), N = 1000 (A) N = 2000 (o), N = 3000 (□), N = 5000 (o). 

[4] Spin Susceptibility Xs for both regular (o) and fiuctuating (o) lattices vs number of 
nodes N. 

[5] Topological Susceptibility Xv for both regular (o) and fiuctuating (o) lattice vs 
number of nodes N. 

[6] Power (o) vs exponential (□) fits for a regular lattice N = 2888 at /3 = 1.7. 

[7] Spin Correlators for regular and fluctuating lattices at criticality together with 
power law and exponential fits respectively. 

[8] Power (o) vs exponential (□) fits for a dynamical lattice N — 3000 aX (3 — 2.3. 

[9] log (xs) vs t~ together with KT linear fit. 

[10] Fluctuation in local scalar curvature and cross correlator of vorticity with loop 
length W vs P . 

[II] Histograms of joint probability distributions for vortex- vortex (Pyy (r)) and vortex- 
antivortex {Pyav (r)) pairing for (3 = 3.0 and /? = 1.0. 
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Topological Susceptibility and Specific Heat 

N=968, regular lattice 
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